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Abstract. We present results for the lowest-order non-vanishing quark number susceptibilities. These were calculated using 
an improved action viz. the HISQ action, which controls taste violations that are responsible for distorting the light meson 
spectrum. Our calculations, with a pion mass of 160 MeV, are also much closer to the physical limit than previous studies. 
We find a broad crossover from the hadronic to the quark regimes, although interactions remain significant even at the 
highest temperatures studied. Our results are also in good agreement with Hadron Resonance Gas models below the crossover 
temperature. 
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The phase diagram of QCD at non-zero values of the quark chemical potentials {jig > 0) is a subject of intense 
investigation. Progress is hampered by the fact that ab initio numerical simulations (i.e. the lattice) are cuiTently not 
possible in this domain. On the other hand, from a variety of model studies it is believed that a second-order phase 
transition exists at some value of the baryochemical potential jig. This picture is currently being tested in the low- 
energy runs at the Relativistic Heavy Ion Collider (RHIC) at Brookhaven [1,2]. 

The basic observables in these experiments are the baryon and charge distributions, their fluctuations and higher 
moments. At vanishing chemical potential these moments may be calculated on the lattice by the method of Taylor 
expansions. These moments provide valuable information about the nature of the QCD chiral transition (which is a 
crossover for physical values of the light quark masses) as well as about the freezeout conditions. Furthermore, in the 
vicinity of the critical point these moments are expected to exhibit scaling behavior. Given sufficient statistics, such 
behavior is in principle observable and can serve as a signal for the elusive critical point. 

SECOND ORDER SUSCEPTIBILITIES 

The susceptibilities X2 ™d x'u (where i,j — u,d,s) are defined as 



r2 d{^i/TY 



and 



r- d{ni/T)d{n,/T) 



(1) 



P and T are the pressure and temperature, while the /x,'s are the quark chemical potentials. The diagonal susceptibil- 
ities measure fluctuations in the quark number distributions while the off-diagonal susceptibilities xW measure 
the correlations between different flavors. By a straightforward change of basis, it is possible to obtain from these 
susceptibilities the corresponding ones in the baryon number {B), electric charge (Q) and strangeness {S) basis. These 
latter susceptibilities are the ones measured in experiments. 

Our calculations were performed with the staggered action. As is well-known, the staggered action describes four 
equivalent species (called tastes) of Dirac fermions in the continuum and thus the staggered spectrum contains sixteen 
pions instead of one. Taste symmetry is broken by interactions at finite lattice spacing and as a result fifteen of the pions 
become heavier by 6{asa^\ This distortion of the light spectrum has consequences for thermodynamics, for e.g. the 
transition temperature is shifted upward. It also affects the calculation of the susceptibilities since these are sensitive 
to the light degrees of freedom. Fortunately a great deal of progress has been made toward designing actions that 
minimize these effects at a given lattice spacing. Currently the state-of-the-art in this respect is the HISQ action [3], 
which is also the action we used for our calculations. 



We calculated these susceptibilities on 2 + 1 -flavor configurations at 10-15 temperatures in the range 125 < T < 340 
MeV. The temperature was varied by changing the coupling and at each coupling the quark masses were also tuned 
so as to stay on the line of constant physics i.e. the mass of the lightest Goldstone mode (identified with the pion) 
was fixed at approximately 160 MeV, while the mass of the corresponding strange meson (identified with the kaon) 
was fixed at the physical value of the kaon mass. The number of configurations used at each temperature varied from 
500-700. 

Our lattice sizes were 6 x 24^, 8 x 32^ and 12 x 48^. At any given temperature, a larger lattice implies a smaller 
spacing (which can be seen for e.g. via T^^ = aNx). These numbers may be compared, for instance, with an earlier 
study by the RBC-Bielefeld Collaboration [4]. That study was done using lattices of size 4 x 16^ and 6 x 24^ i.e. 
farther from the continuum limit. Also, the calculation was done with 2+ 1 -flavor ensembles generated with the p4 
action, which has greater taste symmetry breaking. Lastly, while the kaon mass was fixed at its physical value the pion 
was quite heavy, niji ss 220 MeV. 
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FIGURE 1. Results for X2 ™d X2' compared with earlier results obtained using the p4 action [4]. The black lines marked "SB" 
denote the limiting value of an ideal gas of massless fermions.The transition temperature is ~ 190 MeV for the p4 action at Nx = 6, 
and fa 160 MeV at yV^ = 8 for the HISQ action. 

Fig. 1 shows our results for the lowest-order susceptibilities X2 x|. In comparing these with the earlier p4 results, 
two differences are immediately evident. On the one hand the crossover is smoother in the HISQ case, implying that 
the crossover from hadrons to quarks occurs over a broader range than was previously believed (This is also seen in the 
behavior of the off-diagonal susceptibilities Xbs and Xqs)- Furthermore the approach to the Stefan-Boltzmann limit is 
also slower and both the light and strange susceptibilities never attain this value even at the highest temperatures that 
we studied. 



FLUCTUATIONS OF CONSERVED CHARGES AND HRG 
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FIGURE 2. Baryon Number (X2' ^^f'-) ^^'^ electric charge (X2' right) fluctuations. "/^ scale" refers to the zero-temperature 
observable used to determine the temperature. The black curves are predictions from HRG models. All figures are preliminary [10]. 



The fluctuations of relevance to experiment are those of the baryon number, strangeness and electric charge. If we 
know aU the {u,d, s) fluctuations at a given order, we can compute the {B, Q, S) fluctuations up to that order. Our results 




FIGURE 3. Off-diagonal correlations, between B and Q (left) and between B and S (right). The black curves are again predictions 
from HRG models, and "/k" again refers to the observable used to determine the temperature. All figures are preliminary [10]. 



for the diagonal B and S fluctuations are shown in Fig. 2 while those for the off-diagonal susceptibilities are shown in 
Fig. 3. 

The solid lines in both plots are Hadron Resonance Gas (HRG) curves. HRG models have been found to provide a 
satisfactory description of nuclear matter at temperatures below the crossover [5, 6, 7, 8, 9]. The partition function is 
the product of those for ideal Bose (Fermi) gases, one for each meson (baryon) viz. 

lnir= £ In ^i%{mi,ii,)+ ^ ln^%{mi,^,). (2) 

Jl,K,T],... p,n,A,... 

The low- and high-temperature limits of these susceptibilities have an especially simple interpretation. At very high 
temperatures the degrees of freedom will be quarks. Of these the u, d and s all carry baryon number but only the s 
caiTies strangeness as well and therefore Xn ^ 1 /3- Similarly, at low temperatures the lightest degrees of freedom 
carrying B are the protons and neutrons, but only the protons carry both B and Q. One therefore expects xfp/ ^ 1 /2 
as r ^ 0. 

To summarize, our calculations of the lowest quark and conserved charge susceptibilities have brought forth some 
interesting features. The transition from hadronic to QGP regimes was found to be broader than previously expected, 
and this was also reflected in coiTelations between the conserved charges. The Stefan-Boltzmann limit was not attained 
for T ^ ITc- On the hadronic side good agreement with HRG models was found with the latter providing both a 
qualitative understanding as well as a quantitative description of the data. In the future we shall look to calculate the 
higher-order susceptibilities. These are expected to exhibit scaling behavior and thus should be more sensitive to the 
nature of the phase transition. 
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